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PACS 05 . 40 . Ca - Noise 

PACS 05.45.Xt - Synchronization; coupled oscillators 

Abstract - Common noise acting on a population of identical oscillators can synchronize them. 
We develop a description of this process which is not limited to the states close to synchrony, but 
provides a global picture of the evolution of the ensembles. The theory is based on the Watanabe- 
Strogatz transformation, allowing us to obtain closed stochastic equations for the global variables. 
We show that on the initial stage, the order parameter grows linearly in time, while on the late 
stages the convergence to synchrony is exponentially fast. Furthermore, we extend the theory to 
nonidentical ensembles with the Lorentzian distribution of natural frequencies and determine the 
stationary values of the order parameter in dependence on driving noise and mismatch. 



Introduction. — Synchronization of oscillations by a 
periodic forcing is a general phenomenon observed in nu- 
merous experiments. In this setup the system follows the 
driving and has, in particular, the same frequency, so one 
often speaks on frequency entrainment. Much more non- 
trivial is the effect of synchronization by an external noise. 
Here one can also distinguish between the cases when the 
driven system is entrained by the noise (synchrony) and 
the situations when the noise is not followed (asynchrony) . 
While the difference between these regimes can be hardly 
seen by observing just one responding oscillator, it be- 
comes evident if an ensemble of identical systems driven 
by the same noise is observed: in the case of synchroniza- 
tion all the oscillators in the ensemble follow the forcing 
and their states thus coincide, while in the asynchronous 
state the states of systems remain different. This effect 
is therefore called synchronization by common noise [IHS] ■ 
An interesting realization of this type of synchronization 
is the effect of reliability of neurons [7] . Here one does not 
use an ensemble of identical neurons, but takes one neu- 
ron and applies the same pre-recorded noise to it several 
times. The synchronous case appears as a reliable respond 
to the forcing where all the noise-induced spikes are at 
the same position at all runs, while for asynchrony (an- 
tircliability) the same noise produces different spike pat- 
terns [8][9] . Synchronization by common noise was also ob- 
served in physical experiments with phase-locked loop |10j 



and noise-driven lasers 1 1 1 j . 

Synchronization by common noise can be characterized 
by the largest Lyapunov exponent of the noise-driven dy- 
namics. This exponent governs the growth/decay of small 
perturbations of a synchronous state; a negative exponent 
corresponds to synchrony while a positive one to asyn- 
chrony [IHS] (notice that here the largest Lyapunov expo- 
nent can be interpreted as a "transverse" one, determining 
the growth/decay of the difference between oscillators in 
the ensemble). For periodic oscillators, which in the noise- 
free case have a zero maximal Lyapunov exponent, small 
noise generally leads to a negative exponent (while large 
noise can desynchronize) ; for chaotic systems with a posi- 
tive Lyapunov exponent, strong noise can synchronize (see 
examples of the synchronization-desynchronization transi- 
tion in |l|2|H2Hl4j ) . Calculation of the Lyapunov exponent 
is a relatively easy numerical task, and in many cases it 
can be obtained analytically [JHB]. This theory is, how- 
ever, restricted to the linear analysis of a stability of the 
synchronous state, and does not allow one to follow the 
evolution starting from a broad distribution of the phases. 

The goal of this letter is to present a global analytic the- 
ory of the synchronization by common noise, i.e. a theory 
describing the evolution toward synchrony of the popu- 
lation starting from the distributed, asynchronous state. 
Our theory is based on the Watanabe-Strogatz ansatz [15] 
that allows one to describe a population of phase oscil- 
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lators under common forcing via closed equations for the 
macroscopic, global variables. We will show that the re- 
sulting equations can be written as a noise-driven Hamil- 
tonian system, and will analyze the evolution of ensem- 
bles close to synchrony (where the results of the Lyapunov 
analysis will be recovered) as well as the evolution starting 
from the maximally asynchronous state. 

Global variables description of ensembles. Our 

goal is to describe an ensemble of identical (phase) oscilla- 
tors, and one can show that the following model is general 
enough to describe all interesting cases: 

<p k = n(t)+hn(F(t)e- i '?*), k = l,...,N. (1) 

Here oscillators arc described by their phases ipk, an d Q(t) 
and F(t) are timc-dcpcndcnt common forces (without loss 
of generality we can assume that F(t) is real, otherwise 
one shifts the phases ip by the argument of complex F and 
correspondingly redefines f2). We want to characterize the 
evolution of the ensemble especially for noisy forces. 

Now we present two particularly important applica- 
tions, which can be described using Eqs. (TTJ). One is an 
ensemble of self-sustained oscillators under common exter- 
nal force. Unforced oscillators are described by ipk = oj, 
and a noisy forcing is typically described by 



tpk = lo - cr£(i) sin<p fe 



(2) 



This system, previously considered in [HIS], reduces to 
Eq. (fTJ) with constant f2 = uj and F(t) = cr£(i). Another 
relevant physical setup is a sequential array of shunted 
Josephson junctions subject to a common current I(t): 



2eR~dt 



I c sin tp k = I(t) . 



(3) 



Here ipk is the junction's phase (difference of the phases of 
the macroscopic wave functions in superconductors con- 
stituting the junction), R is the resistance of the shunt, 
I c is the critical current. Supposing that the current I(t) 
has a constant and noisy component I = Iq + 7i(f),thcn 
by rescaling time t —> q2 * R t, we can write system © as 



•4>k 



er£(t) - sin ip k 



(4) 



which is also Eq. (fTJ) with constant F = 1 and time- 
dependent Q(t) = (I + /i(i))(Ic) _1 = w + er£(t). 

In the seminal work [15] Watanabe and Strogatz (WS) 
demonstrated that the ensemble (fT]), for any N > 3, can 
be fully described with three global variables and N — 3 
constants of motion. We will use here the formulation 
of the WS theory given in [IB]. The transformation to 
the global variables p, <I>, \& and constants ipk is performed 
according to 



_ i<i> 



«W>k-*) + 1 



pe 



(5) 



with an additional condition J^k e = 0- The closed 
system of equations for p, $ reads (as 'J does not enter in 



the dynamical equations for p and $, it does not to be 
taken into account.) 



^=L_^ Re (F(t)e- < *), 

$ = m + l±^Im(F(t)e-^) . 
2p 



(6) 



The physical meaning of the global variables p, $ is clear 
from their definition ([5]). The case of uniformly spread 
constants of motion ipk is easiest to interpret, because as 
it has been shown in |16U17| , for a uniform distribution of 
constants ipk one has pexp(i$) = A^ 1 J2k ex P(^fe)- This 
means that z = pe 1 * is the complex Kuramoto order pa- 
rameter widely used for characterizing synchrony in the 
ensemble. Thus, p is roughly proportional to the mean 
field amplitude: for p = the phases ifk are uniformly 
spread while for p = 1 they form a cluster (from which 
at most one oscillator with tpk — & = 7r may deviate) of 
perfect synchrony. The variable $, being the phase of the 
mean field, characterizes the position of the maximum in 
the distribution of phases. Finally, the variable $ charac- 
terizes the offset of the phases of individual oscillators with 
respect to For a non- uniform distribution of the con- 
stants ipk, the complex variable z does not coincide with 
the Kuramoto order parameter; nevertheless, the limits 
p — > and p — > 1 correspond to fully asynchronous and 
fully synchronous cases, so p yields a proper characteriza- 
tion of synchrony. 

Hamiltonian formulation. — Our goal in this paper 
is to describe the evolution of the ensemble (fT]) for noisy 
fi(t), F(t), by virtue of the global variables dynamics ©. 
Remarkably, one can reformulate the basic equations ([6]) 
as a Hamiltonian system |18j . Indeed, in variables 



pcos$ 

"=VT=7 

the equations read 



psin4> 

Vw 1 



(7) 



• dH r./^ i + g 2 + 2p 2 



(8) 



2v/l+p2 + 9 2 ' 



with Hamiltonian 



H(q,p,t)=n(t) 



p 2 +q 2 



F{t) pV±±Z±Z. (9) 



One can also formulate the dynamics in "action-angle" 
variables [15] . where the angle is the WS variable $ and 
the action is defined according to 



J = 



2(1 -P 2 ) 



(10) 
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The Hamiltonian equations now read 
v/2J(2J 
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J = F(ty 



- cos (*) = -^ 



$ = n(t) - F(t) 

with Hamiltonian, 



4J + 1 . ... <9# 

sm($) = —— , 

2y/2J(2J+l) 9J 



(11) 



H(J,<p,t) = Sl(t)J-F(t) 



^2J(2J+1) 



sin$ 



(12) 



The action J yields, according to its relation to the Ku- 
ramoto order parameter p, a natural characterization of 
synchrony in the population of oscillators: J — > cor- 
responds to a maximally asynchronous, uniformly spread 
state, while J — > oo corresponds to a perfect synchrony 
where all oscillators cluster in a single point in phase space. 

The Hamiltonian formulation allows us to give a gen- 
eral qualitative description of the dynamics. Notice that 
since the Hamiltonian is time-dependent, the energy in 
eqs. (|8|11[) is not conserved. Typically, noise leads to a 
growth of energy, either diffusive or exponential (in ex- 
ceptional cases, e.g. if in (fT2l F(t) vanishes, the system 
possesses an integral and no growth of energy is observed) . 
Thus, the action variable J grows and the system tends 
to synchrony. We analyze below this growth for large and 
small values of J, i.e. close to synchrony and close to 
asynchrony, using (J) as the order parameter characteriz- 
ing the level of synchrony. 

Dynamics close to synchrony. Close to syn- 
chrony, i.e. for J> 1, we can approximate the Hamilto- 
nian (32]) as H(J, 0, t) = - Fit) sin $) J which leads 
to a skew system where the dynamics of $ does not depend 
on J: 



j = F(i)Jcos($) , 
$ = Q(f) -F(i)sin($) 



(13) 
(14) 



This yields In J{t) = In J(0)+ f Q F(t') cos $(t')dt' and $(f) 
is a solution of (fT4|) . On average, In J grows linearly in 
time with the rate given by the Lyapunov exponent 



dt 



In J) = -A = (F(f)cos($)) 



(15) 



The same Lyapunov exponent appears when one directly 
analyses stability of the cluster ip\ = ip% = . . . = ip^ = $ 
in system ([1]). The phase of the cluster obeys (fT4|) and the 
small deviation of one of the phases 5<p obeys 

-j^Sip = —F(t) cos($)<5( y 9 , 

what means that Stp ~ exp(At) with the Lyapunov expo- 
nent defined in (fl~5j) . For both applications ([2]) and (j4|), the 
calculations of the Lyapunov exponent for a white Gaus- 
sian noise £ (t) have been already reported in the literature 
(see Fig. 9.4 in book for (3} and refs. gHE| for ©)• 
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Fig. 1: (Color online) Evolution of {J(t)} for noise-driven os- 
cillators for different a (right panel: the same in a rescaled 
time), for ui = 2 and different noise intensities a. Dashed line 
corresponds to relation (11911 . 



Dynamics close to asynchrony. Here we de- 
scribe the dynamics for small values of J, i.e. close to 
the asynchronous regime with nearly uniform distribution 
of the phases pk ■ We perform the analysis separately for 
the noise-driven oscillators ([2]) and noise-driven Josephson 
junctions ([!]). 

An ensemble of noise-driven oscillators ([2]) is described 
by Hamiltonian © with F(t) = a£,(t) and fi = uj = const. 
For small g,|)C 1 we can approximate the Hamiltonian 

as, 

H(q,p,t)=cj^±^ + aml- (16) 

Assuming for simplicity of presentation that we start from 
the vanishing order parameter, i.e. q(0) = p(0) = 0, we 
can easily solve the resulting linear equations: 



? (f) + ip(t) = - / exp[iw(t' - m(t')dt' 



(17) 



which yields for the action J, after averaging, 

(J{t)) = \(q 2 +P 2 ) = 

= ^-f Ioe i ^ t '- t "Hamt"))dt'dt" 

Introducing r = t' — t" and assuming that (£(f')£(^")) 
depends only on r one can integrate over t' so that 

(J(t)) = fy £ (e(O)e(r)) cos W r(l - ^)dr . (18) 

Asymptotically, for large t this describes a diffusive linear 
growth of J (provided the integral converges, i.e. the cor- 
relation function of noise decays fast enough) ; in the case 
of white noise, when (£(0)£(f')) = S(t'), we get 



(J(t)) = -s-i 



(19) 



In Fig. Q] we compare for model ([2]) the theoretical results 
with the numerical ones (obtained from eqs. . (|TT|) ) . 

In the case of noise-driven Josephson junctions ((4]) we 
need to analyze the Hamiltonian 



7J(J,$,t) = (w + a£(t))J- 



y / 2J(2J+l) 



sin $ 



(20) 
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Fig. 2: (Color online) Evolution of {H — Ho) for the ensemble 
of Josephson junctions (j4| for to = 5 and different a (right 
panel: in a rescaled time). The dashed line corresponds to 
formula lt22l). 



We restrict our attention to the nontrivial case uj > 1, 
so that the Hamiltonian is bounded from below. The 
"lowest-energy" state here for a = is the steady state 
with $ = tt/2, J = (u - Voj 2 - l)/(4\/u; 2 - 1). This 
steady state does not correspond to full asynchrony (which 
is not a stationary solution of the equations, because the 
phases rotate non- uniformly) , but to a stationary distri- 
bution of the phases of the ensemble (j4]). Close to this 
equilibrium, we can linearize the equations of motion, so 
that the Hamiltonian in the vicinity of Jo, $o reads 



H{dJ,6<f>,t) = H +a£{t)SJ ■ 



. (21) 



with wi = 0.5^2J (2J + 1), w 2 = 2(uj 2 - 1)(2J (2J + 
l))" 1 / 2 . The solution of the linear equations of motion 
(starting from the equilibrium point) is, similar to (|17[) . 



Ul 1 ? 

<S$ + i— 6 J = a 

K 



exp[in(t' - t)]£(t')dt' 



where k 



1. After averaging we obtain 



(8$f + ^{5J) 2 ) = 



(mC(r)) cos kt(1 - Ll)dT . 



This relation means that, asymptotically, the "energy" de- 
fined as H — Ho — 0.5(wi ((5$) 2 + w 2 (<5J) 2 ) grows linearly 
in time, for the white noise we get 



(H - H ) 



U)l<T 



t 



(22) 



This relation is checked in in Fig. [5J where numerical sim- 
ulations of eqs. (f5|). (fTTj) are presented. 

The results above show that there are two stages for 
the transition to synchrony in ensembles of oscillators 
driven by common noise. As the system is described by a 
noise-driven Hamiltonian, it is natural to characterize the 
evolution through the "energy" - the value of the noise- 
independent part of the Hamilton function. At the initial 
stage, close to asynchrony, the growth of the energy is dif- 
fusive, its averaged square grows linearly in time according 



to (fTi?)) and (|2"2")h When the energy reaches a level of order 
one, a crossover to the other type of behavior, namely to 
an exponential growth of energy, is observed from numer- 
ical simulation (cf. Figs. 1 1 121) . This latter stage means 
that the final convergence of the ensemble of oscillators to 
a synchronous cluster is exponentially fast. 

Nonidentical oscillations. — Here we extend the 
theory to the case of nonidentical oscillators ([2]) having 
a distribution of frequencies g(uj). In this case one first 
generalizes the WS description by assuming the frequency 
dependence of the variables p(u>),^(u>) [TBJ. It is conve- 
nient to introduce one complex variable z(oj) = pe 1 * that 
obeys 



(23) 



The distribution of the phases is now characterized by 
the global order parameter Z = J z{u)g{uj)dijj. Following 
the approach of Ott and Antonsen |20j . it is possible to 
obtain a closed equation for this order parameter in the 
case of a Lorentzian distribution g(uj) = tt 7(7 + (oj — 
wo) 2 ) -1 - Then, assuming that z(cj) as function of complex 
uj does not have singularities in the upper half-plane, one 
can perform the integration to get Z = z{ujq + i"f). Thus, 
the equation for Z follows from (|23[) with u — > luq + 27: 



Z = iujoZ — 7Z - 



(24) 



Transforming from Z to the canonical variables as in 
(|7|) . (|TU)) . we obtain the same equations as © and ([TJ"T) 
but with additional non-Hamiltonian, damping terms: 

l + q 2 + 2p 2 2 2 

q = uj oP + t7g(t) r — - 79 1 + V + 9 ) , 

2 a/1 +p 2 + q 2 

p = -c q ^(t)^—— - 7P(1 + P 2 + I 2 ) , 

(25) 



and 



j = a £( t ) ^ 2J( y + ^cos($) - 2 7 J(l + 2 J) 



$ = ojq - cr£(i) 



4J+ l 



2^27(27+1 



sin($) 



(26) 



With damping terms, the energy does not grow indefi- 
nitely, but saturates as shown in Fig. G3 The saturation 
level corresponds to a bunch of oscillators that do not form 
a perfect cluster, but have a finite spread. 

The stationary level of the order parameter can be esti- 
mated for a state close to synchrony and for a wide distri- 
bution (asynchrony). Close to synchrony, i.e. for J 1, 
we have instead of (|T5|) 



dt 



In J) = cr£(t)cos$ - 27 J 



(27) 



If we neglect the fluctuations of the growth rate and as- 
sume <t£(£)cos<I> = (cr£(t) cos $) = —A, then the station- 
ary value of J is J st = | A|/(2-y). Close to asynchrony we 
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Fig. 3: (Color online) Evolution of the order parameter accord- 
ing to Eqs. (P5|) . ptjl) for oj = 2, 7 = 1CT 3 and different noise 
intensities a. 

can use the same approximation as in (fT6|) , but now the 
equations read 

q = oj Q p - jq + 0.5a£(t) , p = -uj Q q - -yp , 

with the average stationary energy 

a 2 f°° 

which in the case of the white noise yields (J) = a 2 / (I67). 

Remarkably, in both limits the average value of J scales 
as 7 _1 . This is confirmed by numerics presented in Fig. 2) 
In this figure we present also simulations of the oscillator 
populations, which fit nicely the results from the modeling 
of the WS variables. 

Discussion. — In this letter we have developed a 
global theory of synchronization of oscillator populations 
by common noise. Our analysis is based on the Watanabe- 
Strogatz ansatz |15j . which is not restricted by a number 
of elements in the populations and results in an explicit 
time-dependence of the global variables on the common 
forcing terms. These variables can be interpreted as order 
parameters characterizing the population of identical os- 
cillators. For noisy forcing we thus obtained a closed set 
of stochastic differential equations for the global variables. 
An important step in our consideration is a representation 
of the WS equations as a nonautonomous Hamiltonian sys- 
tem; transition to synchrony then appears as the growth of 
energy due to the noisy driving. While for the situations 
close to synchrony the results are essentially the same as 
the previous ones derived from the linear perturbation ap- 
proach, we have demonstrated that when starting from a 
broad initial distribution, the energy first grows linearly, 
and only after a formation of a concentrated cluster the 
exponential convergence to synchrony sets on. Further- 
more, by virtue of the Ott-Antonsen theory [5D] we have 
extended the analysis to populations of non-identical os- 
cillators with a Lorentzian distribution of natural frcquen- 



Fig. 4: (Color online) The stationary level of the order pa- 
rameter for an ensemble of the oscillators with a Lorentzian 
distribution of frequencies, vs. the distribution width. Param- 
eters: u)q = 2, a = 0.3. Squares: modeling eqs. (|25|l . (|26|l : filled 
circles: modeling the ensemble of N = 5000 oscillators for the 
same parameters. 

cies. Here the theory is valid in the thermodynamic limit 
of very large ensembles only. 
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